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ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ 


УДК 621.95.08:51-74 
В.Л. ЗАКОВОРОТНЫЙ, ФАМ ДИНЬ ТУНГ, НГУЕН СУАН ТЬЕМ 


МОДЕЛИРОВАНИЕ ДЕФОРМАЦИОННЫХ СМЕЩЕНИЙ ИНСТРУМЕНТА 
ОТНОСИТЕЛЬНО ЗАГОТОВКИ ПРИ ТОЧЕНИИ 


Дано математическое моделирование упругих деформационных смещений вершины режущего инструмента и 
обрабатываемой заготовки в точке их контакта. Приведены алгоритмы и результаты идентификации пара- 
метров модели для случая продольного точения. 

Ключевые слова: динамическая система, металлорежущий станок, процесс точения. 


Введение. Исследованию и моделированию деформационных свойств подсистем режущего инст- 
румента и обрабатываемой заготовки уделяется неизменное внимание [1-4]. Это связано с тем, 
что при анализе устойчивости процесса резания и автоколебаний необходимо, прежде всего, 
иметь модель деформаций вершины режущего инструмента относительно заготовки. Аналогичная 
проблема стоит и при изучении точности обработки, особенно в случаях, когда заготовка имеет 
значительные деформационные смещения, изменяющиеся вдоль траектории движения инстру- 
мента относительно заготовки. Во всех случаях деформационные смещения как вершины режуще- 
го инструмента, так и заготовки в точке контакта с ней вершины режущего инструмента опреде- 
ляются по отношению к несущей системе станка. В рассматриваемом случае — по отношению к 
его направляющим. Сложность вычисления деформационных смещений вершины режущего инст- 
румента относительно заготовки заключается в том, что они формируются в результате накопле- 
ния деформаций в пространстве всех конструктивных элементов подсистем инструмента и заго- 
товки, расположенных между несущей системой станка и рассматриваемыми точками. При этом 
большое влияние на деформационные смещения оказывают свойства сопряжения конструктивных 
элементов. Конструктивная сложность и неопределенность деформационных свойств узлов со- 
пряжения приводит к тому, что при математическом описании динамики процесса резания необ- 
ходимо не только знать дифференциальное уравнение динамики, но и разработать методы иден- 
тификации всех параметров этой модели. 

В общем случае для анализа динамики процесса резания используются пространственные 
конечномерные модели, приводящие к необходимости анализа следующего дифференциального 
уравнения [4]: 


а4?Х аХ 
О, ТОН =Р(Х,5,1,)+ Л@, (1) 


где т(Х,5,,,)= м, ,(Х,5, 6), ИСХ, = АС, 5.) (5) = 4,8) 


соответственно функциональные матрицы инерционных и диссипативных коэффициентов, 
а также функциональная матрица формирования упругой составляющей сил в зависимости 
от вектора деформационных смещений и технологических режимов. Размерность матриц 


5,К=1,2,..,6; Х=Х,Х,,Х.,Х.Х.,Х} - вектор упругих деформационных сме- 


щений вершины режущего инструмента (первые три координаты) и заготовки в точке 
контакта с ней режущего — инструмента (последние три координаты); 


Е={Е(Х,5, 1), Е(Х,5 >), ВХ, 1), (Х,5 1), 5 (Х,5 >>), В (Х,5, и)" — век 


тор-функции динамической характеристики процесса резания, раскрывающие зависимость 
сил резания от упругих деформационных смещений инструмента и заготовки, а также от 
технологических режимов: величины подачи на оборот 5, и глубины резания 1, при за- 


данной скорости; (г) = {1 (1), Л, (#), №, (#), №. (0, №. (0, (0 = изменяющиеся во времени со- 
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ставляющие сил резания, не объяснимые в координатах упругих деформационных смеще- 

ний, которые интерпретируются как шум. 
Нами поставлена задача математического моделирования, разработки алгоритмов иден- 
тификации и на этой основе — исследования свойств только одной функциональной матрицы 


сх, 5 и.) = [6.4 (Х,5,ы,)] . 


Постановка задачи. На основе результатов исследований, а также с учетом анализа работ, 
проводимых многими учеными, можно сформулировать следующие наиболее важные свойства 
деформационных смещений для подсистемы режущего инструмента. 

1. Все модели должны учитывать пространственные деформации. Представление зависи- 
мости деформационных смещений от сил в скалярном виде, например, в виде одного параметра 
технологической жесткости, является неправомерным, так как в этом случае, в частности при 
изучении механизмов потери устойчивости, теряется информация о циркуляционных силах, прин- 
ципиально влияющих на потерю устойчивости стационарных траекторий движения инструмента 
относительно заготовки. 

2. Эксперименты показывают, что при приложении к вершине инструмента внешней по- 
стоянной силы неизменного в пространстве направления существуют такие направления ориен- 
тации силы, которым соответствуют деформационные смещения этого же направления. Такие 
деформационные смещения являются коллинеарными. В общем же случае, кроме коллинеарных 
деформационных смещений, существуют ортогональные смещения, нормальные к коллинеарным. 
На рис.1 приведен пример диаграммы экспериментально полученных деформационных смещений 
в плоскости, нормальной к оси вращения заготовки для станка 1К62. 
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Рис. 1. Типичная диаграмма суммарных деформационных смещений вершины режущего инструмента 
в зависимости от направления и модуля внешней силы 


На приведенной диаграмме стрелками показаны направления смещения ортогональных 
составляющих деформационных смещений (штрихпунктирная линия) по отношению к их колли- 
неарным составляющим (пунктирная линия). Из диаграммы видно, что существуют два направле- 
ния приложения силы в рассматриваемой плоскости, которым отвечают только коллинеарные на- 
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правления. Это направления, обозначенные «К». Угол между этими направлениями всегда равен 
п/2. Точечными кривыми показаны деформационные смещения, соответственно равные силам 
100 и 200 кг. Важно подчеркнуть, что при одном и том же знаке силы только направления кол- 
линеарных смещений (направления «К») практически не меняются. На приведенной диаграмме 
координата Х’, соответствует направлению, ортогональному оси вращения заготовки, координата 


Х, соответствует направлению скорости резания. 


3. При изменении знака внешней силы, особенно в области малых деформационных сме- 
щений, становятся заметными гистерезисные явления. Если ориентация сил меняется в пределах 


от —20° до +150°, а также при условии, что деформационные смещения рассматриваются в ва- 
риациях относительно некоторого установившегося предварительно напряженного состояния, то 
гистерезисные смещения практически отсутствуют. Эти случаи соответствуют естественным на- 
правлениям изменения сил резания в вариациях относительно точки равновесия, т.е. случаям, 
которые рассматриваются при изучении динамики процесса резания. 

4. По мере увеличения сил неизменного направления наблюдается тенденция уменьше- 
ния приращения деформационных смещений, вызванных одинаковыми приращениями сил. Таким 
образом, для малых деформационных смещений в окрестности равновесия их матрицы жесткости 
можно считать постоянными. Кроме этого, указанные матрицы являются симметричными, т.е. по- 
тенциальными. 

Приведенные свойства позволяют рассматривать деформационные смещения подсистемы 
инструмента в окрестности равновесия как идеальное упругое тело, а зависимость между силами 
и деформационными смещениями моделировать в виде: 

сх =Е, (2) 


ГД а ПЕРО е=[е, + | 5,К =1. 2,3. Причем матрица с= [| — ве- 


щественна в силу гипотезы об идеальности упругого тела, симметрична и положительно оп- 
ределенна, т.е. потенциальна. 


Для оценивания элементов матрицы с = Е, | в качестве Исходной информации принима- 


ют траектории деформационных смещений вершины инструмента относительно несущей системы 
станка в зависимости от сил, имеющих в пространстве неизменную ориентацию. При этом исполь- 
зуются алгоритмы скользящей линеаризации, рассматривающие оценки 


с(х°) = ВЕ, ОХ, жихак = а” $ =1,2,3. (3) 


Придавая силам измеримые приращения ЛР“ ={АА,0,0}', ЛЕ“ ={0,АЕ,,0}", 


при:Х=Х” 


АЕ = {0,0, АЕ} "и измеряя вариации И” деформационных смещений, получаем системы для 
оценивания элементов матрицы с(Х") 
х°с°(Х“)=АЕ°, (4) 
и», "0: 0“ 0. 30-0. -0 
0-30. 0 оса 0 10. 3 
о бое 
о о д 10 0.0 -9 
где = 0 0 0 п хх” хм х 0 0 01| - матрица экспериментально опре- 
оо бо бо 
а > 9 9..0 50 00. 0 
0. 9 бе о и. © 9 
20 40;> 505 027 20. 2 д" 


деленных деформационных смещений; лЕ” = {дЕ®,0,0,0, ДЕ?) ,0,0,0, АЕЗ)}" — вектор изме- 


хе хо 
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ренных внешних значений сил, действующих на режущий — инструмент; 
С? = {С1>б>1,631>61,2>С5.25бз2›С1з»С2 3,633} - Вектор подлежащих оцениванию параметров мат- 
рицы жесткости. 
Для фиксированного значения точки равновесия из выражения (4) получаем массив зна- 
чений с, Ь1=Ъ2,3 с,; =с,; При 1* 1, по которому получаем оценки параметров жесткости: 


их математические ожидания и дисперсии. Математические ожидания Е, являются оценками 


матрицы жесткости с в системе координат ОХ,Х,Х, (рис.2) для заданной точки равновесия. 


Для них также определяются дисперсии 6 оцениваемые среднеквадратичными отклонениями. 
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Рис.2. Система координат, в которой отсчитывается упругое 
деформационное смещение и внешние силы 


Моделирование упругих деформационных смещений в линеаризованном виде, т.е. в виде 
матрицы жесткости, позволяет представить деформационные смещения в новой системе коорди- 
нат Оу у›у, (см. рис.2), которая может быть получена на основе вращения системы координат 


ОХ, Х,Х, с помощью углов Эйлера. В новой системе Оу,у›,у, матрица жесткости является диа- 


гональной. 
Для определения матрицы линейного преобразования от системы координат ОХ, Х,Х, к 


системе координат Оу 57: выполним три последовательных положительных вращения на углы 


Эйлера ф, 0, зу вокруг осей ОХ,, Оу, Оу,. В результате получим матрицу преобразования в 


виде: 
с0$ фсозли — зш фсоз05шлу —с05 фзш ли — зш фсо$ 0 созу — зшфзт 0 
А =| $т фсоблу + с0$ фсо$0зшли —зт фз у + с0$фсо$ 9 созлу —с0$ фзш 0 |. (5) 
эт 9 зшлу эт 9 созли с059 


Для обратного преобразования от системы координат Оу уу, в систему координат 
ОХ, Х,Х, можно использовать следующую матрицу: 
с0$ <фсозли — зт фсо5 Озшлу — зш фсо5ли + с0$ фсо$05шлу  зшлузт 0 
Т = | - созфзш лу — зт фс0$ 9 озу —зт фзш у + с0$ фсо5 9 созлу зтб0созлуи , Т=Л". (6) 


зш 9 5шф — с05 фзш 9 со5 9 
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Отметим некоторые свойства матрицы преобразования Л: 1) лЛ'=Л'; 2) Л"л=Е - 
единичная матрица, поэтому матрица Л является ортогональной. Кроме того, скалярное произ- 


ведение векторов (Л’ 5 )=1, 1=12,3, № - вектор, элементами которого являются элементы 


1-го столбца матрицы Л. 
Связь старых координат с новыми определяется формулой: 


Х=Ау, АР=ЛАЕ" или у=ТХ, ДЕ' =ТАЕ. (7) 
Подставляя формулу (7) в выражение (3), получаем: 
су РЕ ЛЕХ с = Ал. (8) 


Если для малых деформационных смещений в вариациях относительно точки равновесия 
оси (7,,7›,/.)’ можно считать главными, то сила, имеющая направление, совпадающее с одной 
из осей (у, у,, у)“, вызывает упругое деформационное смещение только в этом направлении. 
В данном случае матрица упругих деформационных смещений с’ является диагональной. Таким 
образом, задача заключается в том, что необходимо выбрать углы Эйлера в уравнении (8) так, 
чтобы матрица с’ была диагональной. Это можно сделать на основе решения уравнений 


с} (ф,9лу) = 0,1 = 7, вытекающих из (8). Однако эта процедура для пространственного случая 


представляет большие сложности, так как полученные недиагональные элементы матрицы с” из 
уравнений (8) являются сложными функциями от синусов и косинусов углов Эйлера. 

Матрица жесткости с в статике является вещественной, положительной и симметричной. 
Для любой вещественной симметричной матрицы с существует такая ортогональная матрица Л, 
что Л'сЛ =с®, где с” = 4а2(\.,^.,^-) — диагональная матрица, диагональными элементами 


которой являются собственные числа матрицы с, повторяющиеся согласно их кратности [5, 6]. 
Таким образом, можно привести симметричную матрицу жесткости с к диагональной 
форме, и вычислить матрицу преобразования Л на основе определения ее собственных чисел и 


соответствующих им собственных векторов. При этом матрица жесткости с” = Фаз( 1, А.Л) в 


новой системе координат (у,,у›, уз)” имеет диагональные элементы, являющиеся собственными 


числами матрицы с. 

Приведем более важный случай для матрицы жесткости с, когда она имеет различные 
собственные числа. Это естественный случай для подсистемы инструмента. Предположим, что 
матрица с имеет различные собственные числа ^., Л,, ^., являющиеся решением характери- 


стического уравнения 
ае(е-ЛЕ)=0. (9) 
Для каждого собственного числа /., определим собственный вектор <’ = (&,1,8,2, 8,3)", 
который удовлетворяет системе уравнений 
(с-^\Е) 8" =0. (10) 
Собственные векторы 2', 2”, <*, соответствующие различным собственным числам ^., 
^>, Аз, взаимно ортогональны. 
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Найдем по векторам ©', =”, е*° нормированные собственные векторы е',е?,е? при ус- 


ловии (е’,е’)=1, #=1,2,3: 
1 
В, (11) 
1 
Ё 
Тогда ортогональная матрица преобразования Л, столбцами которой являются нормиро- 
ванные собственные векторы е', 








поз 
А=(е,е ,е’)=|е,2 © 62|. (12) 
@з @23 633 


Теперь можно определить углы ориентаций Эйлера путем сравнения элементов матриц 
(5) и (12) 


9 = агссозе, ,, 


с0509 =е. ,, 
' р : 31 
упфят0 =ел, или + ф=агези(-), (13) 
: аи. 1- (езз) 
эт фэшлу =е 3, 
13 


ди = атс). 
усе, 


Деформационные свойства подсистемы режущего инструмента в плоскости. Эллипсоид 
жесткости для пространственного случая преобразуется в эллипс жесткости в плоскости (рис.З3). 
Для этого случая можно сразу получить выражение для ориентации эллипса жесткости и значе- 
ния его диагональных элементов. 





Рис.3. Система координат, в которой отсчитывается упругое деформационное смещение 
вершины режущего инструмента и внешние силы 


Если идентифицированы матрицы упругих деформационных смещений с в пространстве 
Х , то их выражения с“^ в пространстве Л получаем на основе следующего преобразования: 
И (14) 
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Пусть в окрестности точки равновесия Х” ={Х',Х.}" идентифицированы матрицы де- 


формационных смещений х = {х, ‚х.}’. Тогда для этой матрицы выполним преобразование (14) 


поворота координат 
(^)_| 608% —ф5та | си с. | соза эта 
С == ) я = 
эта  с05$а |с. с››|-зшоа с0$9 
р : (15) 
(с11С0$7@ + С > ЗП” а -— 2, › $ 9.608 0) [(с>2 -с1) 9 29+ с. > ©0829] 
= й а 2 р ы 
[(с›2 — с) 29+ с, > 0829] (сз ” а + > > 608 4+ 2с, > 1 0,С0$ 0) 


Определим условия, при которых недиагональные элементы в матрице (15) были бы рав- 
ны нулю [(с. › — с.) т 29, + с, > с0$29] = 0, тогда получим: 





о = —агсю - | (16) 


Матрица с“^? является диагональной, т.е. 


0 
с — Сол (17) 


7 
0 Со 
2 я . . > ы д. 2 р 
ГДе су/ = (с, 1с0$^@ + с > $” 9 — 26, > 1 060$ 0); Со = (С 91” а +, ›60$ + 26, › $1 9, 0$ 9). 


Выполняя обратное преобразование координат {,,у,} в х={х, ,х,}’, получаем вы- 
ражение для определения матрицы деформационных смещений вершины инструмента при задан- 


ной матрице (17): 
| с0$%  5ш а сх © ре —зт " 
с = = 


—зшо с03а | О су) | зта —соза 


(18) 


р 1 . 
с 1с082 а, + С > 5т? а — (сс!) зщ 2а 
0,1 0,2 260,2 бол 


ИЕ ` » 
— (с, с, )зш 2а ст $т? 9+ С > 608? а 
2 0,2 0,1 0,1 0,2 


По мере увеличения внешних сил изменяются значения матриц жесткости. При этом ори- 
ентация эллипса жесткости практически не меняется. 


На рис.4 показаны зависимости жесткостей от координат {у,,у,}’` е Л. Они хорошо ап- 
проксимируются экспоненциальными функциями, т.е. 








А ПИ 
Сто (У) = со + Со | 1-ехр то ; 
Е | (19) 
В о, 21 
С>о(71) = с50 +650| 1-ехр| - | 
Т 
| у 
где с®,с) - значения жесткостей на начальном этапе приложения внешней нагрузки; 


с®,с09 - значения приращения жесткостей по мере увеличения деформационных смеще- 


ний; Т,?,Т,?) — параметры, характеризующие скорость увеличения параметров жесткости. 
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Рис.4. Зависимость элементов матриц жесткости деформационных смещений вершины инструмента по направлениям 
ориентации эллипсов жесткости в зависимости от деформационных смещений по этим же направлениям 


Практически полное соответствие экспериментальных и теоретических результатов связа- 
но с тем, что параметры матриц жесткости в направлениях {у,,/,} определены на основе иден- 
тификации. Однако они идентифицированы при анализе деформационных смещений по направ- 
лениям {Ё ,0} и {0,Ё, }, а сами деформационные смещения вычислены по алгоритмам [см. 


формулы (15), (18) и (20)] при ориентации силы под углом 45°. Наибольшая погрешность, кото- 
рая не превышает 10%, находится в области малых внешних сил. При больших силах (превы- 
шающих 500 Н) погрешность составляет 1-2%. Приведем также данные об основных параметрах, 
характеризующих деформационные свойства вершины режущего инструмента для различных 
станков токарной группы (см. таблицу). Поперечное сечение инструмента определяется следую- 
щими размерами: ширина основания - 30 мм, высота - 26 мм, вылет инструмента приведен в таб- 
лице. 


Параметры упругих характеристик подсистем инструмента 
для различных станков 


























Тип станка в .10-2, ним вв .10-2,нум © ‚ град Вылет инструмента /, мм 
1Кб2 2,9 5,4 10 50 
1К6б2 1,7 5,1 18 120 
УТ16ФЗ 43 67 12 50 
16К20 с УЧПУ №С210 3,7 6,5 14 54 
18В340ФЗ0 2,8 4,5 19 54 
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Данные таблицы показывают, что параметры, характеризующие деформационные смеще- 
ния, зависят не только от типа станка, но и от параметров режущего инструмента, в частности, от 
его вылета из зажимного приспособления. Нетрудно показать, что при этом увеличение упругих 


деформационных смещений вершины режущего инструмента в направлении Х, возрастают су- 


щественно больше, чем деформации стержня инструмента. Здесь вновь приходится считаться с 
тем, что за счет увеличения вылета инструмента возрастает момент вращения, поворачивающий 
всю суппортную группу, а также перераспределяется нагрузка на отдельные элементы. Ориента- 
ция осей существенно влияет на область устойчивости динамической системы резания, рассмат- 
риваемой в вариациях относительно точки равновесия. 

Особенности деформационных свойств подсистемы заготовки. Главными особенностями 
деформационных смещений заготовки являются следующие: 

— зависимость матрицы жесткости от траектории движения суппорта вдоль оси вращения 
заготовки; 

— практически равные значения жесткостей по направлениям ориентации эллипсов жест- 
кости (рис.5). Поэтому, исходя из уравнений (19), матрица жесткости всегда является диагональ- 
ной для любой ортогональной системы координат, расположенной в плоскости, нормальной к оси 
вращения заготовки. В этом случае для определения упругих деформационных свойств подсисте- 
мы заготовки достаточно знать распределение жесткости вдоль оси вращения заготовки; 

— при рассмотрении проблем устойчивости процесса резания и многообразия, форми- 
руемых в окрестностях стационарных траекторий, необходимо дополнительно анализировать не 
только движения центра масс заготовки, но и крутильные деформационные смещения. В этом 
случае необходимо учитывать крутильную жесткость всей подсистемы привода вращения заго- 
товки. 


Х ‚Х» мкм 


105 



































150 300 450 600 750 мм 


Рис. 5. Зависимость деформационных смещений точки на заготовке вдоль оси ее вращения 
(заготовка закреплена в трехкулачковом патроне и поджата вращающимся задним центром). 
Параметры заготовки: длина 750 мм, диаметр 80 мм, сталь 45 


Выводы. Приведенные математические модели упругих деформационных свойств подсистем ре- 
жущего инструмента и обрабатываемой заготовки позволяют учитывать главные особенности де- 
формаций вершины режущего инструмента и заготовки в точке их контакта. При изучении преоб- 
разования траекторий движения суппорта в траекторию формообразующих движений инструмен- 
та относительно заготовки необходимо учитывать деформации вершины инструмента в простран- 


стве, т.е. Для каждой точки равновесия упругие деформационные смещения инструмента должны 
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рассматриваться в виде квадратной матрицы жесткости деформационных смещений размерности 
3х3. При этом для анализа деформационных смещений подсистемы заготовки достаточно иметь 
информацию о распределении скалярной величины жесткости по длине, дополненную данными о 
крутильной жесткости. Приведенный формализм и методика параметрической идентификации в 
равной мере может быть распространена на другие технологические процессы обработки на ме- 
таллорежущих станках. 

Основные результаты исследования получены при финансовой поддержке РФФИ по проекту 
07-09-90000. 
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